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Abstract 

A model-independent method for the determination of Breit-Wigner resonance parameters is presented. The method 
is based on eliminating the dependence on the choice of channel basis by analyzing the trace of the K and T matrices 
in the coupled-channel formalism, rather than individual matrix elements of the multichannel scattering matrix. 
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1. Introduction 

A general problem in theoretical baryon physics 
is to make a connection between resonances that are 
predicted by various quark models and experiment. 
A reasonable way to proceed is by identifying the 
poles of analytic functions that are able to describe 
simultaneously all experimental data in a multiplic- 
ity of existing channels with theoretically predicted 
resonant states. Therefore, properly and uniquely 
extracting resonance parameters from experiment 
is a task of primary importance. We emphasize the 
problem of uniqueness. The work described here 
is motivated by the need to extract Breit-Wigner 
resonance parameters from multichannel partial- 
wave analyses (PWAs) in a model-independent 
way. Many PWAs of similar experimental data 
produce similar partial waves, while the extracted 
Breit-Wigner resonance parameters are often quite 
different. This fact can easily be seen in the Review 
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of Particle Physics [T] by the Particle Data Group 
(PDG). Each resonance in the Review has been 
parametrized in two ways. First, there are Breit- 
Wigner parameters, i.e. the resonance mass M^, 
decay width F^, and branching ratios Xa (the ratio 
between the partial width into channel a and the 
total width) . Alternatively, T-matrix complex poles 
(Re Wp and — 2Im Wp), as well as complex residues 
(moduli |r|, and phases 6) are given in [T]. Breit- 
Wigner parameters obtained in various partial- wave 
analyses vary quite substantially, partly because 
details of those analyses are different; i.e. the num- 
ber and character of the included channels, different 
parameterization schemes, analyticity constraints 
for scattering amplitudes, the choice of background 
models, and the method of unitarization (if at all) 
of the S matrix. However, it is also the case that 
the methods for extracting resonance parameters 
are different: Argand-plot fits [5], Breit-Wigner fits 
with background |3] , direct fits of analysis parame- 
ters |4l5j . or model specific schemes which extract 
T matrix poles [6l7l8] . 

In this Letter we present a model-independent 
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method for extracting Breit-Wigner parameters 
from any unitary multichannel analysis able to 
provide the full T matrix, using the trace of the 
corresponding K matrix. Since all Breit-Wigner 
parameterizations are equivalent at the energy of 
the /v-matrix pole, the parameters obtained using 
this method should be directly compared to those 
from quark models and lattice QCD. In order to 
connect the results of a model-independent K ma- 
trix extraction with those of a model-dependent 
analysis, e.g. based on the T matrix, we shall keep 
the relations defining multichannel T and K ma- 
trices as general as possible. It turns out that the 
T-matrix trace simplifies the formalism without loss 
of generality, and shows resonant behavior more 
prominently than any T matrix element does. To 
illustrate this, we shall take the T matrix from an 
earlier analysis [7] and recalculate the resonance 
parameters. The T-matrix trace happens to show 
resonant behavior at energies matching those of the 
iC-matrix poles. 



2. Multichannel scattering 

The essence of any multichannel theory is the fact 
that the evolution of a system is no longer described 
by scalars, but by operators acting in an orthonor- 
mal wave-function space, and the transition proba- 
bilities for physical (measurable) processes are given 
by the matrix elements of their representation in the 
chosen basis. Once this basis is specified, the evolu- 
tion of the system is described by solving equations 
which are matrices in the multichannel space, rather 
than scalar equations. 

All equations given here are considered to be ma- 
trix relations, unless matrix indices are explicitly 
stated. The transition probability Pa^b that a two- 
body system from initial channel \a;q) ends up in 
the final two-body (or quasi-two-body) channel |6; q) 
is given by the absolute square of the scattering S''^- 
matrix element PQ_tf, = q|iS''^|a; g) p, where (j des- 
ignates all quantum numbers conserved in the scat- 
tering reaction, and a and b are channels. In the 
case of ttN scattering we have conserved spin, parity, 
and almost conserved isospin (charge symmetry is 
only slightly violated). Conservation of probability 
is ensured if the S matrix (for simplicity, we drop q 
henceforth) is unitary. Therefore, the S matrix can 
be written as S = e^"* , where 5 is some matrix Her- 
mitian in the channel indices. Hermitian matrices 
have real eigenvalues and are diagonalized by uni- 



tary matrices. The 5 matrix is related to a real, di- 
agonal matrix 5d by a unitary transformation 5 = 
U^SdU, where U is a. unitary matrix. The S matrix 
is evidently diagonalized by the same transforma- 
tion, so S' = We^^^'^U. 



The K matrix [HE] is defined as A" i(/-S')/(/-h 
S), where / is the unit matrix. The K matrix can, 
in the eigenstate basis, be written using the diago- 
nal matrix Sd &s K = U'' tanSoU . The K matrix 
is Hermitian because S is unitary, and symmetric 
because of time-reversal invariance, so K is, in fact, 
a real matrix. Thus, J7 is a real orthogonal matrix 
that we henceforth designate as O. 

Every diagonal N x N matrix can be spanned in 
the ortho-normal vector basis {E^, . . . , E'^}: 



E' = 



1 




Vo . . 0/ 



/o . . o\ 



■ 
Vo . 1/ 



so, in our case, we have 



N 



tan Sn = "S^ E^ tan (5' 



i=l 



(1) 



(2) 



where S'' is the ith diagonal element of 6d, also 
known as the eigenphase shift, and N is the number 
of channels. We define the coupling matrices to 
be 

X' = O^E'O, (3) 

and these matrices turn out to be ortho-normal pro- 
jectors: 



Ex' 

i=l 



I, X' = X^S^J, 



(4) 



where Sij is the Kronecker S symbol. 

The trace of a matrix is, by definition, a sum of 
its diagonal elements. A trace has two particularly 
important properties: i) the trace of a product of 
matrices is invariant with respect to cyclic permu- 
tations, Tr [ABC] = Tr [BCA]; and ii) the trace is 
a distributive function with respect to scalars a and 
/?, TT[aA + l3B] ==aTr[^] +/3Tr[B]. 

The orthogonal transformation in definition ^ 
conserves the trace of a matrix, so 



Trx^ 



1. 



(5) 



It follows that the matrices K and T can be written 
as the sums 



2 



N N 

K^J2^^ tan 5^, T^J^ e'^' sin , (6) 

where the connection between the K and T matrices 
is given by the relation 

K = T/{I + iT). (7) 



gives the coupUng strength of the resonance with 
mass and total decay width from channel a 
to channel b. The diagonal element of the matrix 
is the branching ratio x"^ of a given resonance to the 
channel a 

< = xL- (12) 



3. Breit-Wigner parameterization 



4. Extraction procedure 



Elements of tan 5d, as well as the X"' i are func- 
tions of energy or a corresponding kinematical vari- 
able, and their description requires modeling of the 
energy dependence of numerous functions. We sec 
resonances in scattering reactions as real poles of the 
K matrix. The rth element of the diagonal matrix 
tan Sd can be written as [5] 



tan 6'' 



JV2_ 

Mr - W 



tan S 



(8) 



where the selected pole term is parametrized in 
Breit-Wigner form, and it is singled out from other 
contributions, designated collectively as the back- 
ground term at resonance tan (5^ . The Breit-Wigner 
mass {Mr) and total width (F^) parameters arc 
allowed to be functions of the center-of-mass total 
energy W . The reported Breit-Wigner parameters 
M^ and are given by the values of Mr{W) 
and Tr{W) evaluated at an energy equal to the 
corresponding resonance mass M^: 



M^^Mr{M^) 



Tr M. 



(9) 



where we have explicitly written A/,, and F,. from 
Eq. (HI) as functions of energy W. 

The corresponding K and T matrices are given 
by the equations 



F' /? ^ 



(10) 



f;/2 



N 



Mr - W -iV'j2 



where the second term in each equation is the 
coupled-channel background contribution, and FJ, /2 
represents F.r/2 + {Mr — W) tan (5^. When W equals 
the mass of the resonance, FJ, is manifestly equal to 
F^. Although these relations are in general a sum 
over several resonances r, here they are written for 
one resonance for simplicity. 

If there is a pole in the K matrix at some energy 
M^, then the matrix clement x^ab that energy 



The channel dependence of resonance parameters 
can be reduced significantly by using only diagonal 
elements of the T and K matrices. In practice, these 
matrices can be obtained cither by unitary coupled- 
channel partial- wave analyses, or by using partial- 
wave T matrices obtained in diverse single-channel 
PWAs as input to a unitary coupled-channel formal- 
ism, and refitting them to obtain a unitary set of all 
coupled-channel T matrix elements. 

Channel dependence is completely removed from 
the sums 



N 



N 



Tt{K) = ^tan5^ Tr(T) = ^ e'*' sinS^, (13) 

because the traces of the K and T matrices are the 
same as the traces of their similar diagonal part- 
ners tan Sd and e**^" sin(5£), respectively. The same 
is also evident from Eq. ([5]). Consequently. Eqs. (fTO|) 
and pT|) are simplified by taking the traces 



N 



Tr(if) YT^^^ + tan 5^ , 



Mr - W 



(14) 



Tr(T) 



n/2 



N 



Mr~W -iT'j2 



^e'^'siny. (15) 



The last relation, i.e. the T-matrix trace, would be a 
good starting point for model-dependent extraction 
methods. However, instead of putting considerable 
effort into modeling the background and energy- 
and channel-dependent resonance parameters, we 
use the following procedure: 

(i) The parameter extraction procedure starts 
when a full T matrix has been obtained from 
an energy-dependent partial-wave analysis of 
experimental data. 

(ii) Contrary to the usual prescription, where 
Eq. (jlip is used to obtain resonance parame- 
ters from the T matrix in a model-dependent 
way, we use Eq. ^ to obtain the full K matrix 
from the known T matrix. 
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(iii) Poles of Tr K are found to obtain a set of res- 
onance masses M^', ■ ■ ■ M^^, where is the 
number of resonances. 

(iv) Multiplying both sides of Eq. HH) by 
(M^ — W) and setting the energy W to the 
value of the fcth resonance mass {Mj}), the 



corresponding resonance width is isolated: 



^ k 



2 hm 



{M^ - W) Tr(i^) 



(16) 



All other contributions to the K matrix 
trace, i.e. background, other resonances, and 
channel-couplings, are removed in this lim- 
iting process (this relation turns out to be 
similar to Eq.(16) in Ref. [10] for the case of 
the various nN isospin channels), 
(v) The branching ratio of a resonance to a given 
channel can be obtained in similar manner, 
but this time using the diagonal iC-matrix el- 
ement, Kaa from Eq. (fTO|) and definition (fT2|l 



— D lini 



(MP - W) Ka 



, (17) 



where, as before, all undesired contributions 
vanish. 

(vi) Steps (iv) and (v) are then repeated for all 
resonances found in (iii). 

5. Results and discussions 

To illustrate the usefulness of our method, reso- 
nance parameters from a unitary, multi-resonance, 
coupled-channel analysis [7] have been extracted. 
The channels used in the analysis were ttA^, r/iV, and 
an effective two-body channel designated as li^N . 
Extracted parameters are given in Table[TJ The pro- 
posed model gives resonance parameters very close 
to the values obtained by a complicated method of 
diagonalizing the matrix of the generalized Breit- 
Wigner function denominator, with minimal calcu- 
lation. 

We have also compared the iiT-matrix trace to that 
of the T matrix. It can be seen in Fig. [1] that the 
Breit-Wigner resonance positions obtained by look- 
ing for the poles in Tr K (indicated by gray vertical 
lines) directly correspond to the positions of peaks 
in Im(Trr), and of zeros in Re(TrT). The peaks of 
the T-matrix elements corresponding to individual 
channels, however, show a certain deviation from 
that behavior. This suggests that fitting individual 
channels in order to obtain resonance parameters in- 



/x /x 

PDG [1] 



[MeV][MeV] [%] [%] [%] 



511 



35-55/30-55/1-10 
1535±20/150±50 

/ 55-90/3-10/10-20 
Vl650±™/150±^o 



Si 



( unknownA 

1 \^a!2090/NE j 



1543 165 39 54 7 

1553 182 46 50 4 

1680 233 64 16 20 

1652 202 79 13 8 

2054 1926 47 3 50 

1812 405 32 22 46 



CO-70/0/30-40 
1440±5[5/350±100 



p / 10-20/6/40-90 
1^1710±30/100±1{5" 



p ( UNKNOWN 1 
11 I a;2100/NE J 



1482 541 61 39 

1439 437 62 38 

1738 170 44 12 44 

1740 140 28 12 60 

2123 379 3 83 14 

2157 355 16 83 1 



Pl3 



10-20/0/>70 
1720±^!;/100±50 



1776 409 

1720 244 



20 

18 



Dl3 



50-60/0/40-50 
1520±1''/120±;5 



/ 50 
I 17i 



1515 121 

1522 132 



56 

55 



60/0/40-50 \ 
700±50/100±50 J 



1818 126 15 15 



n ( UNKNOWN 1 
l'^ I S32080/NE j 



1817 134 9 
2359 1216 26 

2048 529 17 



14 
6 



Di 



40-50/0/50-60 
1675±l''/150±5j; 



Fl5 



00-70/0/30-40 
1680±;0/i3o±i() 



/ UNKNOWN \ 
1" I a;1990/NE J 



Gl7 



10-20/UNKNOWN 
2190±";/450±100 



Gi 



80 

82 



44 

45 

70 

77 

68 

75 



1674 144 36 64 

1679 152 35 65 



1682 144 67 1 32 

1680 142 67 33 



2139 412 7 3 90 

2262 2036 3 2 95 



1806 286 6 94 

2397 1217 16 84 

2125 381 18 82 



Table 1 

Resonance parameters extracted using the _fC-matrix proce- 
dure given in this paper are listed in bold face. The original 
T matrix was taken from Ref. [7] where the channels used 
were ttN, rjN, and an effective two-body channel n'^N. For 
comparison, Brcit-Wigncr resonance parameters from the 
original reference arc shown below. 

troduces an uncontrolled error, which is avoided if 
the trace of the T matrix is used. 

Unexpectedly, and contrary to previous findings, 
the resonance parameters produced by the /v -matrix 
extraction method presented here, are in accordance 
with values obtained by the original analysis as well 
as with the T-matrix trace. The procedure involves 
no fitting, diagonalizing, nor modeling of the energy 
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Fig. 1. The trace of the T matrix and its contributions for partial waves from Sn to G17. The thick black line represents 
the imaginary part (upper graph) and real part (lower graph) of the trace of the T matrix. The thin line is the ttN clastic 
contribution, the dashed line shows the contribution from r}N, while the dotted line gives the effective-channel contribution 
(unitarity channel). Gray vertical lines are plotted at the trace of if-matrix pole positions. 



dependence of the resonance parameters and back- 
ground. Furthermore, a model-independent proce- 
dure cannot be given with the T-matrix formahsm, 
because background makes a substantial contribu- 
tion to the T matrix, even at an energy equal to the 
resonance mass, M^. The T-matrix background is 
removed at a complex energy equal to the T-matrix 
pole position. This might be the reason why extrac- 
tions of T-matrix poles work much better than T- 
matrix extractions of Breit-Wigner parameters. By 
using the trace of the K matrix, background has 
been completely removed from consideration at the 
resonance energies. 

With regard to the differences between the two ap- 
proaches listed in Tablc[Tl it is rather striking that all 
of them can be explained by arguments presented in 
the original analysis. Since an effective ii'^N channel 
was introduced in [7] to parametrize the first inelas- 
ticity in each partial wave, the parameters of low- 
lying resonances should be much better determined 



than those of heavier ones (especially the third res- 
onances in Sii and Dia). A better quality of pa- 
rameters is also expected for resonances that couple 
more strongly to the measured channels considered 
here. Therefore, A^(1720) P13 and the resonance(s) 
in Gi7 have unrealistic parameters since they are 
completely driven by the effective channel, as can 
be clearly seen from Fig. [1] These problems should 
be removed by the explicit inclusion of additional 
channels in the partial-wave analysis. 

The parameters of the two lowest resonances in 
the 5ii and fis partial waves, as well as those of the 
T'ls, second Pn, and Fi^ resonances, are in rough 
accordance with quark-model expectations for their 
masses and partial widths with the exception 
of the mass of the second TJis, which is predicted 
to be roughly degenerate with the second and 
the Z?i5 resonance. This disagreement could be ex- 
plained by the large coupling of this state to the ef- 
fective channel. The large width and the somewhat 
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larger mass of the first Pn (Roper) resonance ex- 
tracted using the iC-matrix procedure bring these 
parameters closer to those of the class of quark- 
model calculations based on one-gluon exchange po- 
tentials and pair creation for strong decays. 



[10] A. B. Gridncv, I. Horn, W. J. Bri scoe and I. I. 

Strakovsky, arXiv: hep-ph/0 408912| v2 (2004). 
[11] S. Capstick and W. Roberts, Prog, in Part, and 

Nucl. Phys. 45 S241 (2000). 



6. Conclusions 



We have presented a model-independent method 
for resonance parameter extraction using the K- 
matrix formalism. It is shown that real poles of the 
K matrix are related to the resonant behavior of 
the trace of the T matrix. Our resonance parame- 
ter extraction procedure is simple and straightfor- 
ward once the full T matrix is known. Unrealistic ex- 
tracted parameters for some higher mass resonances 
point to the need to include additional channels in 
partial-wave analyses. 

At the energies of the matrix poles, the in- 
fluence of background and channel mixing is elimi- 
nated, so only parameter values obtained at this par- 
ticular energy should be compared directly to the 
predictions of quark model and lattice QCD calcu- 
lations. 

This model-independent procedure cannot be ex- 
tended to the T-matrix formalism because back- 
ground makes a substantial contribution to the T 
matrix, even at the resonance energies M^. This 
might be the reason why methods that extract T- 
matrix poles work much better than those which ex- 
tract Breit-Wigner parameters from the T matrix. 
By using the trace of the K matrix, the background 
has been completely removed from consideration at 
resonance energies. 
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